Conductance Characteristics between a Normal Metal and a Superconductor 

Carrying a Supercurrent 
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The low-temperature conductance (G) characteristics between a normal metal and a clean su- 
perconductor (S) carrying a supercurrent I a parallel to the interface is theoretically investigated. 
Increasing I a causes lowering and broadening of (1) coherence peaks of s-wave S, and d-wave S at 
(100) contact, (2) midgap-states-induced zero- bias conductance peak for d-wave S at (110) con- 
tact, and (3) Andreev-reflection- induced enhancement of G within the gap near the metallic-contact 
limit. Novel features found include a current-induced central peak and a three-humped structure at 
intermediate barrier strength, etc. 
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It is well-known that Andreev reflection plays a funda- 
mental role in understanding the transport properties of 
a normal metal/superconductor junction (NSJ) [1]. From 
the current- voltage {I-V), or the differential conductance 
[G(V) = dI(V)/dV] characteristics of the junction, one 
can learn much information about S, includng its elemen- 
tary excitation spectrum and its order-parameter sym- 
metry, etc. Blonder et al. have developed a general 
theory [2] for studying I-V and G(V) of an NSJ that 
allows a dimensionless barrier-strength parameter z to 
range from metallic contact, z = 0, to the tunneling 
regime, z » 1. There only conventional s-wave sym- 
metry for S was considered. Recently, much attention 
has been paid to the conductance characteristics of d- 
wave, cuprate S in both theory and experiment [3-13]. 
Due solely to the sign change of the d-wave gap-function 
order parameter A(k) on the Fermi surface, a zero-bias 
conductance peak (ZBCP) appears in the tunneling spec- 
trum of an N/(d-wave S) junction with non-(nOm) con- 
tact [4,5,6]. The ZBCP arises from a sizable number of 
midgap states formed at the S side of the N/S interface 
and appears for all z but is narrower and taller for larger 
z. In a large magnetic field, the ZBCP splits into two 
peaks [7,8,10,11]. It is interesting to also study the ef- 
fect of a supercurrent I s in S on G(V). Very recently, 
G(V) for tunneling into a diffusive s-wave superconduct- 
ing wire carrying an I s was measured and compared with 
theory [14]. It was shown that the coherence peaks were 
suppressed and broadened with increasing I s , and the 
effect is the same as that caused by a magnetic field. Be- 
cause the width and thickness of the wire were smaller 
than the superconducting coherence length and penetra- 
tion depth, all variations transverse to the wire could be 
neglected, as was the magnetic field generated by I s . The 
positions of the coherence peaks in G(V) were found to 
practically not shift with I s , up to ~ 4/5 of the critical 
current. In this work, we investigate theoretically the 
conductance characteristics of a clean NSJ with an I s 
in S parallel to the interface by extending the theory of 



Blonder et al. [2]. Contrary to Ref. [14], this work is not 
limited to large z. We consider both s-wave and d- wave 
S with (100) and (110) contacts. Some novel results are 
obtained, especially for z <~ 1, when G(V) does not sim- 
ply reflect the thermally-smeared quasi-particle density 
of states. Hopefully, these predictions can be confirmed 
experimentally. Unlike Ref. [14], the present work does 
not consider Coulomb blockade, which is presumably not 
so important in an extended system and in the clean 
limit. As in Ref. [14], we also assume a uniform I s , and 
neglect self-field. 

When a uniform I s passes through a conventional 
three-dimensional s-wave S, the phase of A(k) has a spa- 
tial variation of 2q s • x, where x is the center-of-mass 
position of a Cooper pair, q s = (m*/2)v s , with v s the 
supercurrent velocity, and m* the mass of a Cooper pair. 
(h = 1 is assumed throughout this work.) At tempera- 
ture T = 0, the magnitude of the order parameter A q 
stays unchanged until the Landau criterion is satisfied 
(i.e. q = 0.5A , where q = q s /k F and A = A /E F ). 
Here Ao is the superconducting gap when I s — 0, kp 
and Ep are the Fermi momentum and energy, respec- 
tively. When q > 0.5A , S becomes gapless, and quasi- 
particles are generated in a portion of the Fermi surface 
[15]. We shall see that this can lead to a ZBCP in G(V) 
for an N/(s-wave S) junction with the barrier strength 
z ~ 1. This current-induced ZBCP is always quite broad 
and not very tall, and its height decreases for larger z. 
It is therefore characteristically different from the ZBCP 
induced by the midgap surface states in a d-wave S with 
non-(nOm) contact that is narrower and taller for larger z 
[4] . (The midgap-states-induced ZBCP has been ubiqui- 
tously observed in high-T c cuprate and other unconven- 
tional Ss.) 

As q is increased further, A g gradually decreases to 
zero at q = 0.67A [Fig. 1(a)]. The supercurrent den- 
sity quickly reaches a peak (the thermodynamic critical 
current density) at q = q c = 0.515A [Figs. 1(c)] [16]. 
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FIG. 1: Dependence of the superconducting order parameter 
on the normalized supercurrent-velocity parameter q for (a) 
an s-wave and (b) a d-wave S. (<f> is the angle between the 
supercurrent and the antinodal direction in the latter case.) 
In (c) and (d), the corresponding dependences of supercurrent 
density on q are given. 



The region q > q Cl in which I s is a decreasing function 
of q, is unstable and can not be observed experimentally. 
(For a two-dimensional s-wave S, superconductivity dis- 
appears immediately after the Landau criterion is met. 
Then q c = 0.5 A .) 

Different from that in an s-wave S, the A q -vs-q rela- 
tion in a d-wave S also depends on the direction of the 
supercurrent. (Here A q denotes the maximum gap in the 
presence of I s .) For a two-dimensional d-wave S with a 
supercurrent, the gap- function order parameter at T = 
is described by [17] 



AO r 

Trln— = J d0cos 2 (29)ln(g + Vp^T), 



(1) 



where g = fL|2g^l|, A 9 = A q /E F , is the angle 
between the supercurrent and the antinodal direction, 
and the integral in Eq. (1) is from to 2ir with the 
constraint g 2 — 1 > 0. 

Figure 1(b) shows the dependence of the d-wave A q 
on q at <f> — and 7r/4. We can see that when q is 
less than ~ 0.3 A , the changes of the order parameter 
with q in both the antinodal and nodal directions are 
almost the same. However, a great difference exists for 
larger q. When I s is applied along the antinodal direc- 
tion, A 9 has a sharp drop (from 0.883A to 0.588A ) 
between q = 0.384A and 0.385A . After that it drops 
continuously to zero at q = 0.53A . When <fi = 7r/4, A q 
gradually decreases to 0.689 A at q = 0.469 A , and has 
no solution beyond. Fig. 1(d) gives the corresponding 
dependences of the supercurrent density on q [17]. It is 
seen that the thermodynamic critical current is reached 
at q = q c = 0.35A (0.39A ) for current in the antinodal 



J 



(nodal) direction. 

The elementary excitations in S are governed by the 
time-independent Bogoliubov-de Gennes equations [18]: 

Eu(x) = h u(x) + J dx'A(s,r)v(x') , (2a) 
Ev(x) = -h v(x) + J dx'A*(s,r>(x') , (26) 

where s = x— x', r = ~(x+x'), and h — ~^1+U5(x)—(J, 
with [i the chemical potential. It is useful to express the 
superconducting order parameter in the form: A(s, r) = 
J dke lk s A(k, r)e l2qs ' r [3]. Neglecting the proximity ef- 
fect near the N/S interface at x = 0, we have A(k, r) = 
A(k)0(x), where Q(x) is a step function, and A(k) is 
the order parameter of a bulk S in the presence of I a . 

In the WKBJ approximation, Eqs. (2) have special 
solutions of the form 



(3) 



where u(x) and v(x) obey the generalized Andreev equa- 
tions [1]: 
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q s • ki 



2m 



)u 



i(k F + q s ) 



q s -k f ^ i(kp-q s ) 



VM + A(k F )9(x)w, 
(4a) 

Vv + A*(k F )e(x)u. 

(46) 

Obviously, the eigenenergy E is symmetric about E = 
q s -kp/m. When q s is applied parallel to the interface 
of the NSJ, i.e. q s = —q s e y , we have 
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(for x > 0) , (5a) 
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(for x < 0) , 



(56) 



where v = sign(kp x ); a v = [—isq 2 /2 + mA^/lkpx], 
with A v = ^/{E + q s k Fy /mf- A J/ (k F )A*(k F ); 
P„ = mi/[-q2/(2m) + E + q a k Fy /m]/\k Fx \; 
7„ = -mis[q 2 s / (2m) + E + q s k Fy Jm]/\k Fx \\ u^ (<) 
and v„^ <S> are constants. For example, in S, we have 
B v = u> /v> = A v {k F )/(E + q s k Fy /m - vA v ). 

Following Ref. [2], after a tedious but straightforward 
calculation, we obtain the Andreev and normal reflection 
coefficients, a(E) and b(E): 



a(E) 



2q + (k + + k-) 



+q+ + 2imU)(k + - g_ + 2imU) - B + (k + + q+ + 2imU)(-k- - q~ + 2imU) 



(6a) 
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b(E) = 



B + (k_ + q_ - 2imU){-k + + q + - 2imU) + B_(k + - q_ + 2imU)(k_ + q + - 2imU) 
5_(-fc_ +q+ + 2imU)(k + - q- + 2imU) - B + (k + + q+ + 2imU)(-k- - q~ + 2imU) 



(66) 
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Here q + = \k Fx \ + + , g_ = \k Fx \ + j + , and k v = 
\kpx\ + volv The critical supercurrent velocity is much 
less than the Fermi velocity. So the Andreev approxima- 



tion, q± « k± « \k Fx \, also holds in the presence of a 
supercurrent. The normalized conductance can then be 
calculated according to a formula given in Ref. [2]: 
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FIG. 2: The normalized differential conductance vs voltage 
for a normal metal/s-wave superconductor junction: (a) z — 
0, (b)z = 0.5, (c) 2 = 1.0, and (d) z = 5.0. Red: q = 0, 
green: q = 0.3A , and blue: q — 0.515A , at which the 
thermodynamic critical current is reached. Note that q = 
q s /k F , and A = A /E F . 



where |Asj? x | — k F cos0, f(E) is the Fermi distribution 
function, G n and G s are the differential conductance for 
S in the normal and superconducting states, respectively. 

S-wave superconductor. In this case, the superconduct- 
ing order parameter A„(kp ) = A q is independent of v. 

In Fig. 2, G(V) at various q and z = 2mU/k F is plot- 
ted. [We have used k B T = 0.01E F and A = 0.1E F .} 
When z = and q = 0, electrons incoming with all mo- 
menta k F with k Fx > can enter S and equal number 
of holes at opposite momenta are retro-reflected into N 
if \eV\ < Ao. So the normalized conductance G = 2.0 
within the superconducting gap if T = 0. With increas- 
ing q, the range of G = 2.0 diminishes and the G(V) 
curve turns into a nearly triangular peak centered at zero 
bias [Fig. 2(a)]. At large z [Fig. 2(d)], the coherence 
peaks are suppressed and broadened with increasing q, 
but contrary to the case of a diffusive superconducting 
wire [14], here the peaks of G(V) move outward while the 



gap shrinks. The intermediate- z results are even richer in 
behavior [Figs. 2(b) and (c)]: A fairly broad and not very 
tall peak appears at zero bias and a three-humped struc- 
ture can also appear for nearly critical q. Note that the 
larger is z, the lower is this current-induced ZBCP. The 
area under this peak is also not conserved as z changes. 
These features are characteristically different from the 
ZBCP induced by the midgap surface states in d-wave S 
with non-(nOm) contacts. [4] 

For electrons entering an NSJ at a fixed incident an- 
gle 8, a ZBCP would result from their contributions to 
the normalized conductance if 2q\ sm9\ > A is satisfied. 
Thus, one can see this peak only if q > 0.5A is satis- 
fied. For 0.5A < q < 0.67A , there is a critical angle 
\6 C \ = arcsin(A°/2g), which decreases from 90° to 48.3° 
in this range. No ZBCP is induced by electrons with in- 
cident angle \6\ < \9 C \. However, only a small portion 
of this regime can be observed, because only the region 
q < 0.515A is stable. 

D-wave superconductor. In this case, the pair potential 
has the form A u (k F ) — A g cos(26' I/ ). Here, = 6 + vex, 
a is the angle between the antinodal direction and the 
positive x axis. 

Figure 3 presents the normalized conductance at dif- 
ferent z and q for a d-wave S with (100) contact (i.e. 
a = 0°). For z = [Fig. 3(a)], the central peak due to 
Andreev reflection is gradually suppressed and slightly 
broadened. For large z [Fig. 3(d)] one sees mainly the 
filling up of the central dip with only a slight suppression 
of the coherence peaks as q increases. For intermediate 
z [Figs. 3(b) and (c)], one sees intricate behavior with 
some similarity to the corresponding figures in Fig. 2. 

Figure 4 shows the normalized conductance at different 
z and q for a d-wave S with (110) contact (i.e. a = 45°). 
It is seen that the ZBCP induced by the midgap surface 
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FIG. 3: The normalized differential conductance vs volt- 
age for a normal metal/d-wave superconductor junction with 
(100) contact: (a) z = 0, (b)z = 0.5, (c) z = 1.0, and (d) 
z = 5.0. Red: q = 0, green: q = 0.2A , and blue: q = 0.35A , 
at which the thermodynamic critical current is reached. 
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FIG. 4: The normalized differential conductance vs volt- 
age for a normal metal/d-wave superconductor junction with 
(110) contact. The z values considered are the same as in Fig. 
3. The q values considered are: Red: q — 0, green: q = 0.2A , 
and blue: q = 0.39 A , at which the critical current is reached. 



states is suppressed, broadened, and eventually split at 
sufficiently large z when q is increased. 

In conclusion, we have studied the differential conduc- 
tance of a clean normal metal/superconductor junction 
carrying a supercurrent parallel to the junction inter- 
face, for barrier strength ranging from metallic-contact 
to the tunneling regime. In the tunneling regime, we 
obtain results similar to the case of a diffusive s-wave 
superconducting wire studied recently, viz., suppres- 
sion and broadening of the coherence peaks for both 
an s-wave superconductor and a d-wave superconduc- 
tor with (100) contact, except that the coherence peaks 
are found to move outward in the s-wave case. For d- 
wave superconductor with (110) contact we also find the 
midgap-surface-states-induced ZBCP to be suppressed 
and broadened and eventually split with increasing su- 
percurrent. In the metallic- contact limit, supercurrent 
causes the Andreev-reflcction-induced conductance en- 



hencement within the (maximum) gap to become weak- 
ened and broadened. For intermediate barrier strengths 
some novel features are revealed including a current- 
induced zero-bias peak and a three-humped structure 
near the thermodynamical critical current density. It is 
hoped that these predictions can be observed experimen- 
tally. We conclude with the remark that this formulation 
can also be applied to the case of an d+s superconductor. 
Because the critical current for an s-wave superocnductor 
is larger than that for a d-wave one, the existence of an s 
component can be verified by a supercurrent reaching a 
magnitude between the critical values of the two waves. 

We wish to thank J. Wei for helpful discussions. This 
work was supported by the Texas Center for Supercon- 
ductivity and Advanced Materials at the University of 
Houston and by the Robert A. Welch Foundation (Ting). 
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